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' Abstract: 

> . 

A method is proposed for the determination of the unitarity angle 
a through tree penguin interference. The modes needed would 
be of the form_B° /B° -> p°M and B°/B° -> uoM where M 

^ ! is spin-0 uu/dd meson, for instance Ai = 7r°, 77, 77', ao or f . 

An analogous method can also determine 7 using M. = K$ or 

^ ■ Kl- The validity of the theoretical approximations used may be 

tested by over determining a with several modes. If two or more 
Qh! modes are used, the determination has a four-fold ambiguity but 

^ e— 1 additional information from pure penguin decays or theoretical 

estimates may be used to reduce the ambiguity to a, a + ir. The 
method as applied to determining 7 is probably less promising. 

h ; 

The early indications of CP violation JI], [J in the neutral B system 
is an important development in our understanding of this phenomenon. It 
is expected that in the near future the angle /3 of the unitarity triangle 
will be determined with considerable accuracy. This is a crucial first step 
in the program of verifying that the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix M of the Standard Model is the origin of CP-violation in the B and 
K systems. To complete this task, the more difficult angles a and 7 also need 
to be determined. Information about these angles, when combined with the 
information concerning the sides of the CKM triangle 0, will either provide 
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an impressive verification for the CKM paradigm or convincing evidence for 
the presence for new physics. 

In this Letter we propose a method to extract a through the interference 
of the b — > uud tree with the b —>■ d penguin; an analogous method may 
also be used for 7 through the interference of the b — > uus tree with the 
b —>■ s penguin. In particular, we will require time dependent observations of 
B° / B — > pM. and B° / B — > ujM. where M. is any self-conjugate spin-0 state 
with a suitable quark content. For instance to obtain a we can use A4 = tt°, 
T], t]', a or /o as we U as related excited states, while to obtain 7 we may 
use M. = K s , K L or related spin kaonic resonances. Comparison of pM. 
and uiM. data for each (pseudo)-scalar M. gives up to a 4-fold degeneracy 
in the determination of sin la (or sin 27) and when information from two 
or more modes is combined, a single value of this quantity should emerge. 
The remaining four fold degeneracy in a (7) can then be reduced to two- 
fold by either using theoretical information concerning the tree or penguin 
amplitudes or by measuring a pure penguin mode related by SU(3). 

There have been several methods proposed to extract these angles. For 
a one can consider oscillation effects in B° — > 7r + 7r~ although one must ac- 
count for the penguin through isospin analysis |J by observing B° — ■> -k°7t°. 
Since the branching ratio to 7r°7r° is expected to be small and hard to ob- 
serve, it may be preferable to consider three tt final states. In J7J, a is 
determined through isospin analysis of B — > tt ± 7i + ti^, B^ — > -k ± -k°'k q and 
— > 7T + 7r~7r°. In this approach one can take advantage of resonance 
effects in the Dalitz plot; however, there may be problems in precise mod- 
eling of the resonance structure. Another method for extracting a from the 
interference of w-penguins with i-penguins in may overcome 

the disadvantages of the 2tt and 3tt final states || although the analogous 
B s decays are required for the analysis. The angle 7 may be extracted via 
direct CP violation through the interference of b — > ucs and b — > cus and 



via time-dependent studies of CP violation in the B s system ]T0 . 

In the method discussed here, we rely on the following two approxima- 
tions: 

(1) The contribution of the electro- weak penguin (EWP) is small. 

(2) The qq pair which arises in a strong penguin does not form the u in 
the final state. 
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Recall that the EWP are also assumed to be small in some of the other 
proposed methods || 0] mentioned in the preceding paragraph for extracting 
a since they rely on isospin. Our second assumption is true at lowest order 



in perturbation theory by color conservation and, as discussed in [11], using 



renormalization group improved perturbation theory one can show that it is 
valid to a few percent in general. 

Let us first consider the case of determining a through the interference 
of the b — > d penguin with the b — > uud tree. For a given final state of this 
form / let us define A to be the amplitude for B° — > / and A the amplitude 
for B — > /. The SM amplitude for B° —+ f can be written as: 



A f = T(f)v u + P u {f)v u + P c {f)v c + P t (f)v t = f f v u + PfV t (1) 

where T(f) is the tree contribution to the amplitude, Pi(f) is the penguin 
contribution due to the diagram with an internal quark of type it, and i>j = 
VitVid- Unitarity of the CKM matrix allows us to express the amplitudes in 
terms of T) = T(f) + P u (f) - P e (f) and P f = P t {f) - P c (f). 

If we denote tf = \TfV u \ and Pf = \PfVt\ then, using the phase convention 



of fll2| , we can write: 



A f = t f e +i ^+p f e- ip 

A f = tfe-'+pfe^ (2) 

In the limit that AT /V is small, which is a good approximation for the B° 
meson, the time dependent decay rate is: 

1 = -e~ |r| (Xt + bY f cos x b r -bZi sin x b r) (3) 

T B o dr 2 \ j j j ) 

where b = +1 (-1) for B° (5°)-mesons, r = Tst and the coefficients Xf, Yf 
and Zj are related to the amplitudes in eq. (Q) by: 

X f = (\A f \ 2 + \A f \ 2 )/2 
Yf = (MM 2 - \A,p/2 

Zj = Im(e-^A}A f ) (4) 
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From these, we can also determine, up to a two fold ambiguity, the quan- 
tity: 

Zf = Re(e- 2if3 A}A f ) (5) 

since {Zj) 2 + {Zf) 2 = X 2 — Y 2 . We can expand X and Z in terms of the 
tree and penguin amplitudes as: 

— \Pf\ 2 + \tf\ 2 ~ 2r/ cosa 
Zf = \pf\ 2 + \tf\ 2 cos 2a — 2rj cos a 
Zj = \tf \ 2 sin 2a — 2rj sin a (6) 

where r/ = Re(tfp*f). 

If we eliminate |tj| 2 and rj from eq. @ and denote = Zf + = 
e~ 2il3 A* f A f , we obtain: 

X f - Zf cos 2a - Z 7 sin2a = X f - Re(Z f e~ 2ia ) = (1 - cos2a)| P/ | 2 (7) 
Likewise, if we eliminate \pf\ 2 and r/ we obtain: 

X/ -Zf = (1 -cos2a)|t / | 2 (8) 

Our primary approach to obtaining a is to use eqn. (0) to relate two 
modes by matching B° — > p°Ai and B° — > uA4, M. = any spin-0 uu/dd 
meson. The fact that the magnitude of the penguin is the same for these 
two modes in principle provides enough equations to cleanly determine a. 
This, however, will normally result in discrete ambiguities, so it is important 
to supplement this method in such a way as to over determine a and thus 
reduce these ambiguities. This may be accomplished in a number of ways: 

(a) Use the same procedure with a number of different candidates for A4, 
indeed M. = 77, 7/, ao, fo} and higher excited states may all be 
viable. In addition a meson with non-zero spin may be used for A4 if 
one analyzes the angular distributions of the final state to determine 
the time dependent magnitude of the specific helicity amplitudes. Of 
particular interest in this category are / = p°p°, p°uj and uouj all three 
of which share a common \p/\ 2 . 
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(b) Use SU(3) to determine the magnitude of the penguin amplitude from 
a pure penguin process. As we shall see, even if there are considerable 
theoretical errors in the application of SU(3), this additional informa- 
tion will often resolve the ambiguity between a and tt — a which cannot 
be resolved purely from the above p — u matching. 

(c) Other theoretical bounds such as the ratio between the tree and the 
penguin, again even with the presence of appreciable uncertainty, may 
be used to constrain the results providing a means to distinguish be- 
tween the different ambiguous solutions. 

(d) Finally, this method may be used in conjunction with other methods 
for obtaining a. In the Standard Model all such methods should agree 
on a single value for a and so all independent methods to determine a 
should be used in parallel. 

We now illustrate our approach to extracting a from a comparison of p°Ai 
with ujM. for the specific case when Ai = 7r°. 

It follows from approximations (1) and (2) then that the penguin con- 
tributes to these final states only via the dddd channel, therefore p p o v o = 
~Vwk° an d thus from eqn. (0): 

X p o n o - X^o = Re(e~ 2ta (Z p o 7V o - Z^o)) (9) 
If we then define 

£ = Xpo^o — X^o; (e td = Z p o n o — Z^o (10) 

we obtain 



a = (6 ± cos" 1 i)/2 + <^ or (11) 

^ [ 7T 

In general, this will give a 16-fold ambiguity in a caused by the four fold 
ambiguity in eqn. ( |TTD together with the sign ambiguities in Z^ and Z^ v0 . 
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However, if |£/C| > 1 for some choices of the signs for Zf, there will be 
no corresponding solutions for a and in such cases the degeneracy will be 
smaller. In addition, notice that if one takes the incorrect sign choice for 
both Zjp^o and Z^ and also takes the incorrect sign choice for cos -1 
one will obtain 7r/2 — a. There will therefore always be a four fold ambiguity 
in this method with a, a + ir, n/2 — a and 3n/2 — a being possible solutions. 
Note that all these solutions have the same value of sin 2a so, in effect, this 
method determines the value of sin2a with up to a 4-fold ambiguity. 

The p — uj comparison is clean in that no theoretical assumptions are 
required beyond (1) and (2) so the determination of sin 2a is limited only 
by statistics. Furthermore it is self-contained in that only information from 
the B° is needed, therefore it may be used at asymmetric B-factories or at 
hadronic B facilities although the ability to detect n° is required in the modes 
we consider. 

As suggested above, the first step to resolve the ambiguity is to apply 
eqn. to as many final states as possible. By combining the results from 
several values for M. , one should be able to extract a unique value for sin 2a 
(hence a four- fold ambiguity for a). 

One way to further reduce this ambiguity to 2-fold is to use SU(3) to 
estimate the magnitude of pj from a related pure penguin b — > s transition. 
In the case of the p°7r° /uix final states, the appropriate pure penguin decay 
is B s — > KgKg or KlK* l where K* SL means a neutral K* meson which 
decays to a final state 71°Ks,l- Note that in this mode the amplitude for 
the spectator to form the pseudo-scalar and the amplitude for the spectator 
to form the vector states are combined in the same way as they are in the 
B° — > p°7T case. This decay is a pure penguin assuming that the rescattering 
contribution of tree diagrams is negligible and so its decay rate allows us to 
estimate \pf\. Generally the penguin amplitude obtained by p— uj comparison 
will be different for a and 7r/2 — a. 

Bearing in mind the SU(3) assumption used to determine \pf\ from B s — > 
K S K* S or K L K* L , this pure-penguin method could also be used as a "stand- 
alone" method for determining a by rearranging eqn. ([?[) into: 



X, - \p,\ 2 = Re((Z f - wv 2 "") 



„ =(w ±co.-'*Ld2£)/2 





or 

7T 
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where Rf and pj are defined by Zf — \pf\ 2 = Rfe lflf . This gives up to an 
8-fold ambiguity: 2-fold from the sign of Z^, two fold from the sign of cos -1 , 
and two fold from the mod n ambiguity. In this context both / = pir and 
/ = uj7c° may be used so the ambiguity of this stand-alone method may, in 
principle, be reduced to the 2-fold mod ir ambiguity. 

Theoretical input can also be used to reduce the ambiguities and refine 
the determination of a. Estimates of the penguin or tree amplitudes may 
serve this purpose. In this context it is useful to note that the ratio of eqn. (0) 
with eqn. @ gives us the tree-penguin ratio as a function of a: 



Pf 



Xf ~ Zf (13) 



Xf — Zj- cos 2a — Zj sin 2a 



A theoretical range for this ratio will therefore translate directly into a range 
of values for a given the experimental inputs for each final state / separately. 

The degree of precision which may be achieved depends very much how 
the different ambiguous solutions happen to align up for various modes and 
in turn on the phases for the various amplitudes involved. In order to illus- 
trate the situation, let us construct a toy model based on related observed 
branching ratios taking M. = tt° and 77. In the case of M. = 77, the analogous 
b — > s penguin modes: B~ — > K*~r\ and B° — > K*°r] have been observed at 
CLEO 0, |l4l with an average branching ratio of ~ 2 x 10~ 5 . We thus ex- 
pect that \p prj \ 2 ~ |^d/Vf S | 2 Br(K*rj), (using amplitudes in units that square 
to branching ratio), so taking |Kd/Vts| ~ 0.2 we obtain \p pv \ 2 = 0.8 x 10~ 6 . 
In order to estimate the penguin rate for p°7r° let us consider the analogous 



pure penguin mode (4>K) which has been observed at BaBaR [181 at a rate 
of Br(B — > <f)K) « 1CT 5 . As before we estimate Ippo^ol 2 ps 0.4 x 10" 6 . 

The tree process in this channel which have been observed at CLEO [19 
are 7i~p° and ir°p~ both of which have a branching ratio ~ 10~ 5 . These are 
color allowed while the processes we are interested in are color suppressed, 
hence we will assume that \tf\ 2 ~ O(10 -6 ). 

For the purpose of generating specific illustrative examples, we will as- 
sume that for each Ai, we will parameterize the amplitudes in the following 
general form: 

t pM =A^ t ulM =a^e^A^ p pM = - PujM = a^e^A* 4 (14) 
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where we will specifically look at the parameter set: 



Al = A% = 10- 3 , al = 1.1 = al = 1.1, a* = 0.65, aP p = 0.90, 
C = -90°, ^ = 15°, ^; = 20°, ^ = 10°, (15) 

which are consistent with the estimated rates above. Further, for the purpose 
of this illustration, we will also assume that the true value of a = 75°. 

Given perfect experimental data, the set of solutions for a in the case of 
M = 7T° are {52.4°, 75°, 75.5°, 114.1°} together with angles related to these 
by a — > 7r/2 — a; a — > 3n/2 — a and a — > 7r + a. In the M. = r\ case, the 
corresponding solution set is {75°, 78.9°, 82.6°, 86.9°} and the related angles. 
Only the true solution and the three other related angles (i.e. 15°, 195° and 
255°) are common to both sets. 

To get an idea of how well we can do with the statistical errors from a 
finite amount of data, let us define: 



No = (number of B°) + (number of B°) ■ (acceptance) 



(16) 



Let us now derive a chi-squared function, % 2 , using statistical errors of the 
input quantities Xf, Yf and Zj. If we define Nf = (Br(B° — > /) +Br(B — > 
f))N /2 then: 



(AXtf = X}/N f (17) 

If we want to determine Yf and Zj from the time dependent distributions in 
eq. (0), we can use the expectation value of time dependent operators which 
are proportional to these quantities. Using the optimal observable as defined 
in |2(J we obtain: 



9y = o ^i+aS* J cos{x b r); g z i = b y I sm^ 6 rj; 

< g Y >= Y f /X f - < g ZI >= Zj/X f (18) 

Assuming a tagging efficiency T, the statistical errors in these observables 
are: 
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< (A^) 2 > 
< (Ag ZI ) 2 > 

In Fig. 1 we show the minimum \ 2 as a function of a for these two cases 
given Nq = 10 9 and tagging efficiency T = 0.5 assuming that the central 
values are those given by the parameters in eq. ([T5|). The M. = tt° case is 
indicated by the dashed curve, the M. = r\ case by the dotted curve and 
the sum by the solid curve. By construction, in each of the cases the curve 
must hit for possible solutions. This is clearly true in the case of A4 = tt° 
while for M. = 77, the solutions occur in tight clumps so, for instance, in the 
range 75° — 86.9°, x 2 is close to but rises rapidly outside of this range. 
Looking at the sum, la uncertainty range about 75° is 75tl\t (as well as 
the corresponding range around 15°). It is asymmetric due to the M. = rj 
curve. If we lower Nq to 2 x 10 8 (i.e. scale \ 2 down by a factor of 5), then 
the range becomes 75lg 5 g 6 where most of the bound now depends on the 
M. = r] case. Likewise for Nq = 10 8 the range is 75lg 8 g 9 . Clearly if a single 
mode has a tight clump of solutions as exemplified by the A4 = r] case, early 
data can put relatively tight bounds on the solution but it is difficult to 
distinguish between the members of the clump. To refine the solution with 
more statistics a mode with more widely scattered solutions such as M. = tt° 
is helpful. 

Fig. 1 also shows how a pure penguin mode may be helpful in distinguish- 
ing between the solution at 75° and the solution at 15°. The magnitude of 
the penguin amplitude in the M. = vr° case for the minimum \ 2 solution at 
each a is shown with the dash-dotted curve. This is clearly very different in 
the two regions illustrating how pure penguin data from B s — > K$Kg may 
solve this ambiguity. Of course the entire graph repeats itself in the range 
180° — 360° leaving a mod tt ambiguity. 

Fig. 2 plots the same quantities as Fig. 1 where we have changed ^ to 
190° to illustrate a somewhat different behavior for the Ai = r] case. In this 
case, the solution set is {75°, 80.1°, 101.3°, 106.2°} and we have somewhat 
better determination of the true solution since the M. = r] mode now only 
has one false solution near to the true one. The two false solutions for M. = i] 







Y 2 




\l + 2xl 


X 2 






z] 


TN, 1 


I 2x1 


X 2 



(19) 
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near ~ 105° are eliminated by the Ai = 7r° data although there is a local 
minimum of x 2 i n that vicinity since the latter case has a false solution at 
114.1°. In particular, for N = 10 9 , the l-o" range about the true solution is 
75^25 while if Nq = 2 x 10 8 then the range is 75^48 but the other minimum 
near 108° begins to become a viable solution. 

The information from this method can and should be combined with 
various other methods for determining a. In particular the method of [|7| 
uses one of the same modes (p°vr°) so jointly fitting a, t p7T o and p p o n o between 
the two data sets will produce more constrained results. In addition, since 
that method uses cross-channel interference, the ambiguity between a and 
ot + ir may be resolved. The interference between i? ± — > p ± p° and B — > p ± u 
in the uj — > 7r + 7r~ channel as discussed in |H| may also be used to address 
this ambiguity. 

The comparison between a as determined here and the value determined 
by the method of |J is particularly interesting since in that case T is purely 
the u-penguin; deviation between the two values of a could therefore be an 
indication of physics beyond the Standard Model. This relation between the 
u-penguin phase and the tree phase may also be tested in the method of [|7| 
through the comparison of charged and neutral B decays. Furthermore, the 
methods in ]lTJ depend only on direct CP violation and so a discrepancy 



with the methods involving time-dependent CP violation could indicate new 
physics in BB oscillation. 

To obtain the angle 7 we now consider modes which are sensitive to the 
interference of the b — > s penguins with the b suu tree. The spin-0 particle 
which recoils against the p/uj should therefore be a K$, Kl or any other spin- 



kaonic resonance that is self conjugate by decaying to a Kg or Kl |§2 
in the final state. For this analysis we need to assume that an accurate 
knowledge of (3 is available, which will likely be the case. 

Let us denote such final states as g. In this case, the weak phase of the 
tree is still +7 but the b — > s penguin phase is ~ in the Standard Model. 
The decomposition analogous to eqn. (0) thus gives: 

Ag = tg e+l1 + Pg 

A g = t g e- l ^+p g (20) 

where i g = \T g V* b V u d\ and p g = \P g V^ h Vtd\- Using the observables previously 
introduced in eqns. (f|,[5]) we can define 
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Z g = Z g e+* f> ; Zf = Re{Z g )- and Z\ = Im(Z g ). (21) 

where we need to know (3 in order to obtain Z from Z. As in the case of a 
we find that: 



X g - Re(Z g e +2 ^) = (1 - cos 7 )N 2 (22) 

We can use this equation in the same way as eqn. (|7|) above since here too the 
penguin produces p and u via the dd channel so we can equate \p g \ 2 between 
ujK s and p°K s . If we now define: 



(e 



ie 



Z P°K S 



-z 



uiK a 



(23) 



then 



^ = -{6 ±cos~ 1 i)/2 or 7i- (^icos- 1 ^ (24) 

In this case, however we are somewhat restricted in the number of modes 
g which are experimentally accessible since spin-0 kaonic resonances are re- 
quired and, aside from K°, these particles are close in mass to higher spin 
resonances with similar decay modes. One could, however, gain additional 
modes by time dependent angular analysis of g = p/u K*° and treating each 
of the three helicity states as a separate mode. 

As in the a case, one can also resolve the ambiguities by comparing 
this to a pure penguin amplitude. In this case the pure penguin which is 
appropriate is — * K°p ± . This has the advantage that it will also be 
produced at asymmetric B-factories. 

Using these methods for the determination of 7, however is likely to be 
less promising than the a determination. Although the branching ratio to 
such modes is relatively large (O(10~ 5 )), the color suppressed tree amplitude 
is only about 4% of the penguin so that the interferences effects required to 
solve for 7 will be < 4%. 
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To see how the method works in this case, it is useful to consider the limit 
that t g « p g , we can then use the following approximate relations: 



tan ^ = I. |2 Z f y +°(\ t 9 /P 9 \ 2 ) (25) 
\Pg I ^9 

tan 7 = - Z Y L ~f K +0(\t g / P gf) (26) 

In the first relation, if \p g \ is known from the analogous pure penguin mode, 
this gives 7 up to a 4-fold ambiguity (2- fold for Z# — > — Zr and 2- fold for 

7 — > 7 + 7r). 

In order to estimate the number of 5-mesons required we note that the 
key measurement would be the determination of Zi which would in turn rely 
on an accurate determination of Zj. Assuming that Zj/X « 0.04, we see 
from eqn. flT|) that to get a 3-sigma determination of this quantity requires 
N g » (X s /A^) 2 (l + 4x 2 )/(2Ta; 2 ) so that if we take Xg/AZ* = 3/0.04 we 
obtain N g = 3 x 10 4 . If the branching ratio for B° — > ujK° is roughly the 
same as the measured branching ratio ||23[ B + — > cjf^ + of 1.5io!e x 10~ 5 , 
then N ~ 2 x 10 9 . With this number of B-mesons, determination of X g 
and \p g \ 2 to a precision of < 1% which is also required should be statistically 
possible. However very tight control over the backgrounds would be necessary 
to accurately determine the denominator of eqn. Q2"5p. Eqn. ( J2"B"D can give 7 
up to an 8-fold ambiguity and the experimental requirements are similar. 

In summary, our primary method for determining a involves observing 
the time dependent decays B°/B — > p°Ai and B°/B — > uj°A4 where M. is 
a uu/dd spin-0 meson. If we assume that electro- weak penguins are negligi- 
ble and that the gluon fragmentation to an u is also negligible, the penguin 
contribution to these two modes has the same magnitude. This allows us to 
solve for a with up to 16-fold discrete ambiguity. These ambiguities may be 
reduced to 4-fold by combining two or more such modes. Using other theoret- 
ical input or estimating the magnitude of the penguin by SU(3) related pure 
penguin modes one can reduce this to a two fold ambiguity. The precision 
that can be achieved depends on the various amplitudes and strong phases 
involved and it was found that reasonable bounds might begin to be achieved 
with iVo ~ 10 8 while with iVo = 10 9 relatively tight bounds on a appear pos- 
sible. In the case of 7 it is possible to use the analog of this method for the 
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decays B°/B — > p°Ks and B°/B — > cji^s and in principle it can be carried 
out with iVo ~ 10 9 . In this case however the relevant interference effects are 
small so that systematic errors in the measurements must be controlled to 
< 1%- 

This research was supported in part by US DOE Contract Nos. DE-FG02- 
94ER40817 (ISU) and DE-AC02-98CH10886 (BNL) 
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Figure 1: The x 2 function for the minimum x 2 solution at various values of 
a is shown for the inputs given in eqn. ( |15|) where the true value of a = 75°. 
The B° — ► p/uj 7T° results are shown as a dashed curve, the B° — > p/uj rj 
results are shown as a dotted curve while the sum is shown as a solid curve. 
The magnitude of the penguin for the minimum x 2 solution is shown in units 
of 10~ 3 by the dot-dashed curve. 
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Figure 2: The x 2 function and penguin amplitude as in Fig. 1 with a = 75° 
and the inputs as in eqn. (^) except with ^ = 190°. 
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